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We study the non-equilibrium dynamics of two component bosonic atoms in a one-dimensional 
optical lattice in the presence of spin-orbit coupling. In the Mott insulating regime, the two- 
component bosonic system at unity filling can be described by the quantum spin XXZ model. The 
atoms are initially prepared in their lower spin states. The system becomes out of equilibrium 
by suddenly introducing spin-orbit coupling to the atoms. The system shows the relaxation and 
non-stationary dynamics, respectively, in the different interaction regimes. We find that the time 
average of magnetization is useful to characterize the many-body dynamics. The effects of even and 
odd numbers of sites are discussed. Our result sheds light on non-equilibrium dynamics due to the 
interplay between spin-orbit coupling and atomic interactions. 

PACS numbers: 03.75.Lm, 03.75.Mn, 05.70.Ln 


I. INTRODUCTION 

Non-equilibrium dynamics of a closed many-body sys¬ 
tem is of fundamental importance [3, i in quantum 
physics and statistical mechanics. For example, the lo¬ 
cal observables appear to be in thermal states even if 
the entire many-body system is in a pure state. The 
eigenstate thermalization hypothesis [2-13 proposed 
to explain such phenomena for complex quantum sys¬ 
tems. However, it cannot be applied to integrable sys¬ 
tems. In fact, a general mechanism of thermalization is 
still lacking |2|. It is important to study the dynamics 
of an isolated many-body system in experiments. This 
may help to understand the behaviour of non-equilibrium 
dynamics and thermalization mechanism. 

Ultracold atoms offer experimental platforms to study 
many-body dynamics of closed systems 0 ■ For instance, 
ultracold atoms can provide a relatively long accessible 
time to study the non-equilibrium dynamics, and the 
high level of controllability to adjust the interaction pa¬ 
rameters for initiating the many-body dynamics. In ad¬ 
dition, ultracold atoms have been exploited to simulate 
a lot of intriguin g q uantum phenomena such as quantum 
phase transition [7| from a Mott insulating regime to su¬ 
perfluid. The advanced detection techniques have been 
invented to enable one to individually address a single 
atom 0 , 11 ]. Therefore, this can be used for probing the 
dynamics in the microscopic description 0- Recently, 
the relaxation dynamics of closed ultracold atomic sys¬ 
tems have been observed PEI- 

In this paper, we consider a system of two-component 
bosonic atoms in a one-dimensional (ID) optical lattice. 
In the Mott insulating regime, the system, with unit fill¬ 
ing, can be described by a quantum spin XXZ model 
[Hi- We consider all atoms to be initially prepared in 
their lower spin states. To study the non-equilibrium 
dynamics, the spin-orbit (SO) coupling [p is suddenly 
turned on. In fact, SO coupling, which produces the 
interaction between the particle’s spin and momentum, 
naturally exists in solid-state materials. It gives rise to 


a number of intriguing effects such as topological insula¬ 
tors and superconductors, etc. PP. Spin-orbit cou¬ 
pling in atomic Bose-Einstein condensates P can be 
produced by inducing two-photon Raman transition us¬ 
ing a pair of lasers [pl- Alternatively, the SO coupling 
between the atoms in the lattice can be induced by pe¬ 
riodically shaking the lattice potential [l^ (pj- More 
recently, the techniques for adjusting SO coupling have 
been shown [U P. The SO coupling gives rise to 
pair interactions between two neighbouring atoms and 
Dzyaloshinskii-Moriya (DM) interactions [23l - l^ in the 
lattice. The DM interaction leads to rich magnetic phase 
diagram [pj, for example, it can induce spin spirals. 

We consider the magnetization as an observable to 
study the many-body dynamics. We find that the time 
average of magnetization is useful for characterizing the 
non-equilibrium dynamics. The system exhibits the re¬ 
laxation and non-stationary dynamics in the different 
interaction regimes which depend on the SO coupling 
strength and the ratio of inter-component interaction to 
intra-component interaction. Indeed, the dynamical be¬ 
haviours relate to the overlap between the initial state 
and eigenstates. The interplay between the SO coupling 
and atomic interactions leads to the changes of this over¬ 
lap, and results in the different dynamical behaviours. 

Thermalization [28i - [^ occurs when a subsystem 
evolves to a mixed state even if the system is in a pure 
state. In the relaxation regime, the spins rapidly relax 
just after the SO coupling is turned on. The degree of 
quantum coherence of local spins can be measured by 
using the purity. An atom in each site will evolve to a 
nearly completely mixed state in a sufficiently long time. 

The transition from the relaxation dynamics to non- 
stationary evolution occurs when the SO coupling 
strength is strong and the inter-component interac¬ 
tion strength becomes sufficiently larger than the intra¬ 
component interaction. In the non-stationary regime, we 
find that the distinct dynamics are displayed for the even 
and odd numbers of sites, respectively. For even-number 
cases, the effective two-level dynamics is shown. This 
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forms a superposition of two distinct states [s^ during 
the time evolution. These superposition states are useful 
for quantum metrology [s^ . On the other hand, the spin 
system becomes ferromagnetic in the odd-number cases. 

This paper is organized as follows: In Sec. II, we 
introduce the system. In Sec. Ill, we study the non¬ 
equilibrium dynamics of this system. We characterize 
the many-body dynamics by using the time average of 
magnetization. We discuss the relaxation dynamics of 
local spins and non-stationary dynamics in the different 
interaction regimes. The effects of even and odd num¬ 
ber of sites on the dynamics are discussed. We provide 
a discussion and a conclusion in the sections IV and V, 
respectively. In Appendix A, we derive the effective Rabi 
frequency. 


II. SYSTEM 


We consider the two-component bosonic atoms to be 
trapped in a one-dimensional optical lattice. We assume 
that this system has open boundary conditions. The two- 
component Bose-Hubbard model can be used to describe 
the interactions of two-component bosons in an optical 
lattice. The Hamiltonian I/bh can be written as, {h = I), 


Hbh = ^ [ja(alai+i -b H.c.) -b ^n“« 

Oi,i 

+Uabnin’l , 


1 ) 

( 1 ) 


where and q;| are the annihilation and creation op¬ 
erators of an atom in the atomic spin state ja), and nf 
is the number operator at site i, and a = a,b. The 
parameter Ja is the tunnel coupling, Ua(b) and Uab are 
the intra- and inter-component interaction strengths of 
atoms, respectively. We assume that the tunnel coupling 
and atom-atom interaction strength of each component 
are nearly equal, i.e., Jq~ J{,~ J and Ua^Ub~U. 

We consider that the atom-atom interaction strengths 
U and Uab are repulsive, and also they are much larger 
than the parameters such as J and tso- In this strongly 
interacting regime with unit filling, it is convenient to 
write the two-mode bosonic operators in terms of angular 
momentum operators, i.e., Sf = a\bi, S~ = b\ai and 
Sf = {b\bi — a\ai)/2. The Hamiltonian can be written in 
terms of spin operators as: 


Hxxz = A 


N-l 

E 

2 = 1 


[ 2(1 - 2Ur)S!S!^,-iS+S-,, + S-Sl,) 


( 2 ) 


where A = 2J‘^/Uab and Ur = Uab/U. The system can 
be described by the quantum XXZ model 321. 

To study the non-equilibrium dynamics, we consider 
the SO coupling to be suddenly applied to the atoms. 
For a non-interacting single-particle Hamiltonian, it is 


given by 


Hqo ~ 2^^ PkxUy -b Skyax, 


( 3 ) 


where k = {kx,ky,kz) is the momentum of a particle 
with a mass m, /3 and 6 are the SO coupling strengths, 
and / and ax^y are the identity and the Pauli spin opera¬ 
tors, respectively. If i 5 eq uals —/3, then it is called as the 
Rashba SO coupling [T^. The Dresselhaus SO coupling 
[13 is an alternative form of spin-orbit coupling where /3 
equals 6 and they are both negative El- Recently, the 
SO coupling with an equal weight of Rashba and Dressel¬ 
haus couplings has been realized in a ®^Rb Bose-Einstein 
condensate [I 4 I, i.e., /3 yb 0 and (5 = 0. By inducing a two- 
photon Raman transition via the laser beams, the equal 
weight of Rashba and Dresselhaus SO couplings can be 
produced, cx kxUy. The Hamiltonian, which describes 
such SO coupling in a ID optical lattice, can be written 

as mill 


Hso = tso - a|5i_i -b H.c.), (4) 

i 

where tso is the strength of spin-orbit coupling. 

In the presence of SO coupling, the effective Hamilto¬ 
nian is written as [ 2 ^, 


N-l 


-i](: 


2Ur-l]SfS^ 




+(^)'( 5 + 5 +1 + s-s-,)-(s+s-^, + S-S+,) 

(5) 


The last terms in Eq. are called the DM interactions 
[ 2 ^ . Now the total system can be described by the 

XYZ spin model with the DM interactions [25l - [^ . 

Let us briefly discuss the various terms in the effective 
Hamiltonian in Eq. ([5]). The first term favors to 

preserve the same polarization of two neighboring spins 
with their initial states. This gives rise to bound magnons 
in the XXZ chain [I^. The second term 
describes the interaction which excites and de-excites the 
two neighboring spins in pair simultaneously. The third 
term + S~S^^ leads to spin-exchange between 

two nearest neighbors. The forth term — S'f 

causes spin-rotation in which the rotation direction de¬ 
pends on the spin states of their nearest neighbors. 

When tso/J and Ur are roughly equal to one, all terms 
equally contribute to the dynamics. All the terms are in 
competition. This results in eigenstates with the different 
combinations of spin states. On the contrary, when both 
parameters tgo/ J and Ur are large, the first term becomes 
dominant. All spins tend to have the same polarization. 
This leads to the very different dynamics in these two 
interaction regimes. 
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III. MANY-BODY DYNAMICS 

We investigate the quantum dynamics of the system 
by suddenly applying the SO coupling to the atoms. Ini¬ 
tially, all atoms are prepared in their lower spin states, 
i.e., 

|vi/(o)) = iu •••;;)• (6) 

This state is indeed an eigenstate in Eq. Q for tgo = 0. 
To initiate non-equilibrium dynamics, a quantum quench 
has to be introduced. By suddenly turning on the SO 
coupling, the system is described by the Hamiltonian in 
Eq. ®. The initial state is no longer an eigenstate of 
the Hamiltonian of the system. If the SO coupling is 
suddenly turned on, then the system becomes out of equi¬ 
librium. We numerically simulate the dynamics of this 
spin chain by using exact diagonalization (see [s^ and 
references therein). 



FIG. 1. (Color online) Contour plot of the average of mag¬ 
netization \M\ versus parameters tso and Ur, where the total 
time taking for average is r = lOA. The system’s sizes A = 11 
and A = 12 are shown in (a) and (b), respectively. 

We study the dynamics of magnetization which is given 
by 

2 ^ 

( 7 ) 

' i=l 

The magnetization M is equal to -1-1 (—1) when all spins 
are in up(down) states. We take the average of magneti¬ 


zation \M\ within a period r = lOA, i.e., 

M=- [ \M(t)\dt. (8) 

'T Jo 

By taking the absolute sign of the magnetization, we can 
ensure that M is positive. Here we set the period r to 
be lOA. This period is sufficiently long until the local 
spins become steady. It enables us to characterize the 
dynamics by using M in Eq. ®. 

In Eig. [l] we plot the time average of magnetization 
M versus the parameters tgo and Ur- When tsoj J ranges 
from 0.75 to 2, M is close to zero. Note that the initial 
magnetization is -1. After taking the average, it becomes 
nearly zero. This means that the system reaches a steady 
state and the magnitude of magnetization is small in a 
long time. 

As both tso and Ur increase, M gradually increases. 
The system is no longer stationary in a long time. This 
shows the different layers in Eig. [T] When both parame¬ 
ters tso/ J and Ur become sufficiently large, M can reach 
nearly 0.5 and 1, respectively, for the even and odd num¬ 
bers of sites. In this regime, this shows that the different 
behaviours of even- and odd-number of sites. 

It should be noted that the average magnetization is 
close to one in Eig. [T]if tso/ J is less than 0.75. In fact, 
we cannot characterize the dynamics of this interaction 
range because the time evolution is too slow and they do 
not reach the steady state within the period r = lOA. 

According to Fig. [U the many-body dynamics can be 
mainly classified into the two different types which are re¬ 
laxation dynamics and non-stationary evolution, respec¬ 
tively. We will discuss them in the following subsections. 

A. Relaxation dynamics 

Now we study the quantum dynamics of local spins 
in the parameter region which the long-time average of 
magnetization M is about zero in Fig. [TJ In Figs. [DJa) 
and (b), we plot the dynamics of magnetization M, for 
the odd and even numbers of sites, respectively. Initially, 
the magnetization M is equal to —I. When the SO cou¬ 
pling is turned on, the magnetization swiftly increases 
in a short time. Then, it becomes saturated to around 
zero in a longer time. This shows that the spins nearly 
relax to the steady states. As the system size increases, 
the magnetization becomes more steady. Besides, the 
relaxation dynamics shows no different in the odd- and 
even-number cases. 

To study the relaxation of local spins, we examine their 
purities. The purity is a quantity which measures the 
degree of quantum coherence of a system. It is defined 
as 

Apur = Tr(p2), (9) 

where p is the density matrix of a system. If the sys¬ 
tem is in a pure state, then the purity is equal to one. 
Otherwise, the purity is less than one. 
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FIG. 2. (Color online) Magnetization M versus time t/X, for 
tao = J and Ur = 1- The odd and even numbers of sites are 
shown in (a) and (b). The different sizes of the system are 
denoted with the different lines: = 11 and 12 (red dotted), 

= 13 and 14 (blue dashed) and A^ = 15 and 16 (black 
solid), respectively. 


We investigate the purity Ap^,. of ion i, where Ap^. = 
Tr(p?) and pi is the reduced density matrix of ion i. The 
reduced density matrix pi can be obtained by tracing out 
the rest of the other spins in the chain. In Fig. [3Ka), we 
plot the purities of local spins versus time. The initial pu¬ 
rity is equal to 1. Afterwards, the purities rapidly drop 
and the purity Ap^. of a local spin decreases to about 
0.5 in a long time. The purity of a spin-half particle in a 
completely mixed state is 0.5, where all diagonal elements 
of Pi are equally weighted and off-diagonal elements are 
zero. This implies that the spins are nearly fully relaxed. 
Additionally, we study the purity Apj)'/ = Tr(p^g^jf) of a 
half of the spin chain in Fig. [3Kb), where the reduced 
density matrix phaif can be obtained by tracing out an¬ 
other half of the chain. The purity becomes saturated to 
a value 0.018. It is very close to the totally mixed state 
which gives the purity = 2 ^/ 2 -Af _ 0.015625 in our 
case. This means that the half of a spin chain can be 
approximately described by a totally mixed state. 

Then, we compare the relaxation dynamics with the 
different strengths of parameters Ur- The average mag¬ 
netization M increases when Ur becomes larger. In 
Fig. SKa), we plot the magnetization M as a function 
of time, for Ur = 1 and 1.5, respectively. The magneti¬ 
zation M in Eq. 0 fluctuates around zero with a larger 
magnitude if Ur becomes larger. Therefore, this will give 
a larger value of M in Eq. ([^ which is obtained by taking 
average of the absolute value of M. We study the purity 



t/x 

FIG. 3. (Golor online) Purities of local spins and half of a 
spin chain versus time t/X. In (a), the purity of local spin i 
versus dimensionless time t/Xis shown. The different lines are 
denoted for the purities of different ion i: i = 4 (solid black), 
i = 6 (blue dashed) and i = 8 (red dotted), respectively. In 
(b), the purity of the left part of a spin chain versus time 
is plotted. The parameters are used: N = 12, tao = J and 
Ur = 1. 


of a local spin in Fig.jTKb). The purity increases when Ur 
increases. This means that this local spin has a higher 
degree of quantum coherence. This suggests that M is an 
useful quantity to characterize the relaxation dynamics 
of this system. 

Indeed, the occurrence of relaxation can be understood 
by examining the overlap between the initial state and 
eigenvectors of the system. We consider the probability 
coeflrcients of the initial state and eigenvectors 

Pn=\{'SmEn)\^ ( 10 ) 

where |'I'(0)) and |if„) are the initial state and the n-th 
eigenvectors. 

In Fig. EJa), we plot the overlap probabilities ver¬ 
sus n, where tso/J = Ur = ^ and n is an index of the 
n-th eigenstate. This corresponds to the previous case 
in Fig. 121 We can see that the initial state has a large 
overlap with the eigenstates. The initial state overlaps 
with almost entire eigen-spectrum. In Fig.j^b), the over¬ 
lap probabilities Pn are plotted versus n. It corresponds 
to the case in Fig. 01 where the magnetization shows 
stronger fluctuations in the dynamics. Obviously, the 
overlap between the initial state and the eigenstates is 
much smaller than that in FigjSJa). Here the off-diagonal 
terms of the observables are suppressed if there is a large 
overlap between the initial states and the eigenstates of 
the system. 
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t/X 

FIG. 4. (Color online) Magnetization M and purity Apu^ 
of the 6-th spin versus time t/X are plotted in (a) and (b), 
respectively, for = 12 and tao = 8J. The different param¬ 
eters Ur ~ 1 and Ur = 1.5 are denoted with the black-solid 
and blue-dotted lines, respectively. 

B. Even-odd effect 

When both parameters tgo/J and Ur are suffi¬ 
ciently large, the dynamics of the system becomes non¬ 
stationary. We find that the dynamical behaviours are 
totally different between the even and odd numbers of 
sites. For even-number cases, the system undergoes an 
effective two-level dynamics. In contrast, the system be¬ 
comes ferromagnetic if the number of sites is odd. 

1. Even number case: Effective two-level dynamics 

When tso/J and Ur are both larger than one, the 
first term in the Hamiltonian in Eq. m becomes domi¬ 
nant, i.e., Hq oc and the other terms in the 

Hamiltonian are perturbations. Obviously, the states 
(I tt • • • tt) ± I II ... ||))/-\/2 are two nearly degenerate 
eigenstates of the Hamiltonian Hq. In this regime, the 
entire many-body dynamics can be effectively described 
by these two degenerate states if the system starts with 
I II ... II) in Eq. (jH]). We plot the time evolution of the 
magnetization in Fig.|6Da). The magnetization shows pe¬ 
riodic oscillations. The effective Rabi frequency decreases 
as the parameter increases and a larger magnitude can 
be attained. In fact, the superposition of the two degen¬ 
erate ground states can be produced, i.e., 

|vI/(t))«.Ci|tt..-tt)+C2|l|...U), (11) 



n 


FIG. 5. (Golor online) Overlap probability versus the 
index n of the n-th eigenstate, for N — 12. In (a), the param¬ 
eters tao = J and Ur = 1 are used. The different parameters 
Ur = I and 2 are used in (b) and (c), respectively, but with 
the same SO coupling, tao = 8J. 


where |cip -I- |c 2 p = 1. 

The two degenerate states can be coupled via the high- 
order virtual transitions. The effective Rabi frequency 
can be approximately obtained which can be derived by 
using the high-order perturbation theory in Appendix A. 
Since this effective Rabi frequency inversely scales with 
the power N, the rate of evolution becomes slow as the 
system’s size increases. This hinders the creation of su¬ 
perposition of two spin states when the system goes large. 

In Fig. [SDc), we plot the overlap between the initial 
state in Eq. ([6]) and the eigenstates. The overlap is much 
smaller than the two previous cases in Figs. EDa) and (b) 
which shows the relaxation dynamics. We have presumed 
that the system is strictly contained in the degenerate 
subspace which can be described by the two degenerate 
states only. Therefore, there are only two degenerate 
states involved in the entire dynamics. In the limit of 
strong interaction, the dynamics cannot be thermalized. 
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FIG. 6. (Color online) Magnetization M versus time t/\ 
for even {N = 12) and odd {N = 11) numbers in (a) and (b), 
respectively, and tao = 8 J. The different lines for the different 
Ur are shown: Ur = '^ (black solid line), Ur = S (blue dash 
line) and Ur = 2 (red dotted line), respectively. 


2. Odd number case: Ferromagnetic 


In Fig.jSKb), we plot the magnetization versus time for 
odd-number cases. The spins tend to maintain in their 
ground states when Ur increases. The spin system is 
ferromagnetic. Indeed, the initial state in Eq. m is an 
eigenstate if the parameters tso and Ur go large. There¬ 
fore, the magnetization is about —1 for the large values of 
tso and Ur- The small fluctuations around -1 are shown 
due to the virtual transitions from the perturbed terms. 

The odd- and even-number cases are totally different 
to each other. In even-number cases, the effective two- 
level dynamics occurs due to the virtual fluctuations of 
the perturbation terms ^ Sj^^^-t-II.c.. However, these 

perturbation terms alter the spin state in pair only and 
therefore they cannot contribute the dynamics between 
the two degenerate states | ft • ■ • t) and | H ... 1) 
for odd-number cases. Also, the DM terms cancel the 
contributions from the spin states and their reflection 
states. For example, the states |001) and its reflection 
states 1100) will be cancelled in the perturbation series. 
Therefore, the DM terms cannot lead to the effective two- 
level dynamics in this case. 


IV. DISCUSSION 


We have investigated the non-equilibrium dynamics of 
a closed quantum system by suddenly applying the SO 
coupling. It is necessary to produce and tune the re¬ 
quired SO couplings to the atoms. There are several ways 
to create SO coupling such as two-photon Raman tran¬ 
sition [13, [II and shaking lattice H^ [2^ . Recently, the 
techniques for tuning SO coupling [2l| have been demon¬ 
strated by using Raman coupling with laser fields. How¬ 
ever, this method may produce unwanted heating due to 
spontaneous emissions of atoms. 

Alternatively, the SO coupling can be exploited by 
shaking the lattice periodically. This “shaking” method 
has been used to successfully generate the artificial gauge 
potential to cold atoms in lattices [s^. More recently, 
the theoretical proposals for the realization of SO cou¬ 
pling have been put forward [H [2^. This method is 
able to create SO coupling without heating if the appro¬ 
priate driving conditions are met M- The other scheme, 
which overcomes the problems of spontaneous emissions, 
has also been proposed [s^. Apart from that, it is also 
required to adjust the inter- and intra-component inter¬ 
action strengths for observing the different kinds of dy¬ 
namics. This can be made by using Feshbach resonance 
[iojl . The scattering length between the different compo¬ 
nents of atoms can be modified by app lying the appro¬ 
priate magnetic helds to the atoms 41|. 

In addition, we make a rough estimation of the re¬ 
laxation time-scale in realistic experiments. We take 
the typical value of the tunnel coupling J to be about 
100^200 Hz in optical-lattice experiments [lOl. l42j| . Since 
the ratio Uab/J is tunable [3], we assume that it ranges 
between 5 to 10 to enter the Mott insulator regime. This 
gives the parameter A = 2J‘^/Uab to be roughly about 40 
to 80 Hz. In Fig. [51 we can see that the system takes 
the period A for relaxation which is about 10 ^ 20 ms. 
To observe the intrinsic effect of relaxation in an isolated 
system, the relaxation time must be much shorter than 
the damping time from the external noise sources. The 
typical heating time of atoms in optical lattice is several 
hundreds of ms 0. This suggests that the relaxation of 
local spins can be detected in experiments. To take the 
time-average of magnetization, the required time is about 
lOA ^ 100 ms which is comparable to the heating time. 
In realistic experiments, the average time of extracting 
the magnetization can be chosen to be shorter than the 
heating time to characterize the dynamical properties of 
the system. 


V. CONCLUSION 

In summary, we have studied the many-body dynam¬ 
ics of two-component bosonic atoms in a ID optical lat¬ 
tice by suddenly introducing the SO coupling. In Mott- 
insulating regime, the system can be described by a quan¬ 
tum spin system. We study the dynamics of magneti- 
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zation of the system. We find that the time average 
of magnetization is useful for characterizing the non¬ 
equilibrium dynamics. The system shows the relaxation 
and non-stationary dynamics in the different interaction 
regimes. In the relaxation regime, the magnetization be¬ 
comes nearly stationary in a long time and the local spins 
become nearly fully relaxed. When the SO coupling is 
strong and the inter-component interaction strength is 
sufficiently larger than the intra-component strength, the 
system becomes non-stationary. The totally different dy¬ 
namical behaviours are shown for the even and odd num¬ 
bers of sites in the non-stationary regime. 


Appendix A: Derivation of effective Rabi frequency 


We study the effective Rabi frequency between the two 
degenerate states from the perturbation theory. We con¬ 
sider the pair-excitation interaction V to be perturba¬ 
tion. The Hamiltonian Hg + V, where Jig and V 

are given by 


N-l 


i7o = A ^ {2[(y)' - l] ( 2 ^ - 


H = A 


i=l 


{-f) {Stst,, + s-s-,,), 


(Al) 


where Hg is treated as an unperturbed Hamiltonian and 
y as a perturbation. 

To calculate the effective Rabi frequency, we need to 
evaluate the virtual transition from | ft ■ • ■ t) and | 

... 4 ,) • The perturbation term V can make the transitions 
for two neighboring spins in pair. It takes N/2 virtual 
transitions from | ft • • • t) and | H ... 1) only. We 
can then make the approximation by using the iV/2-th 
order perturbation theory. However, the DM terms will 
take N virtual transitions to connect these two states. 
The A^-th order perturbation theory has to be used. The 
correction from the DM terms is much smaller than that 
from the term V. Therefore, we can safely ignore the 
DM terms in calculating the perturbation theory if tso/J 
is comparable with Uab/U. 

We can obtain the leading terms of the iV/2-th order 
eigenenergy: 


V i. yi \yi 

j^{N/2) _ ^^2feA7/2-i U li=l ^ ki7 

^ni,n2 ~ / V 

1X2=1 


UZi 

-fother terms, 


(A2) 


where V^l = - E^^ and D de¬ 

notes the degenerate subspace for m = | ft ■ • ■ t) and 
^^2 = I ■ i)) and Nc is the number of possible terms 
that connect the two degenerate states via virtual fluctu¬ 
ations. The number W can be obtained numerically by 
counting all possibilities to connect the two states. The 
leading term of E^^n} can be written as 

~ (A3) 

TABLE I. This table shows the effective Rabi frequencies from 
the numerical results and the perturbation theory for the dif¬ 
ferent system’s sizes, and the percent error r]. 


N 

Q.r/\ 

Dfl/A 

r) 

6 

8.9760 

8.1633 

9.1% 

8 

2.7013 

2.3324 

13.7% 

10 

0.83642 

0.66639 

20.3% 

12 

0.25964 

0.19040 

26.7% 


We compare the effective Rabi frequency between the 
numerics and the approximation from the perturbation 
theory in Table HI for the different sizes N. The percent 
error rj shows the error between the exact numerical value 
and the approximation. It is defined as 


77 = X 100%. (A4) 

Here we denote D/j = E^fn} as the approximation from 
the perturbation theory. This approximation is fairly 
good when N is small. However, as N increases, the 
error grows. In fact, we have taken account of the leading 
term from the A/2-th order perturbation only. When N 
increases, the calculation should include the higher order 
perturbation terms to improve the accuracy. 
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